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The status of the theory of color confinemnt is discussed. 
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1. Introduction 

Quarks and gluons are the constituents of hadrons 
and the fundamental fields of the QCD Lagrangean. 

Lqcd = -\Tr[G^Gn + T, f i>{iD/ - m/)^ (1) 

Quarks can be detected by use of electroweak probes 
at short distances. They have never been observed as 
free particles. The experimental upper limits to the 
existence of free quarks can be found in Particle Data 
Group reports 0. 

The upper limit to the ratio of the quark abundance 
in nature n q to that of the nucleons n p is ^ < 1CP 27 
to be compared with the expected value in the Stan- 
dard Cosmological Model 10~ 12 . 

The inclusive cross section a p to produce a quark 
or an antiquark has upper limit 10~ 40 cm 2 to be com- 
pared with the value expected in the absence of con- 
finement « 10 _25 cm 2 . 

A suppression factor of 10 -15 can only have as a 
natural explanation that n q = and o~ q = due to 
some symmetry. 

This phenomenon is called color confinement : 
asymptotic particle in QCD are only colorless. 

The obvious question is then: does QCD imply 
color confinement, and, if so, by what mechanism? 
What is the symmetry which produces confinement? 

For the first time it was conjectured in ref[2( that 
the Hagedorn limiting temperature Q could in fact 
correspond to a deconfining phase transition from 
hadrons to a gas of quarks and gluons (the so called 
Quark Gluon Plasma). 

Experiments colliding heavy ions have been set up 
at CERN and at Brookhaven to detect such transition. 
It is not clear what would be the smoking-gun signal 
for that transition and no clear statement can be done 
up to now. 

Virtual experiments, i.e. numerical simulations of 
the theory on the lattice have instead demonstrated 
the existence of a deconfining transition. Lattice 
Montecarlo techniques produce a discretized approxi- 
mant to the functional Feynman integral which defines 
QCD . If the lattice spacing a is small compared to 
hadronic scale A and A is in turn much smaller than 
the lattice size La the numerical estimate can be a 
good approximation to the functional integral. 

Also QCD at finite temperature can be simulated 
by similar techniques. The partition function Tr[e~~] 



is equal to the Feynman integral compactified from 
to y in the time direction, with periodic boundary 
conditions in time for bosonic fields, antiperiodic for 
fermions. In formulae 



Tr\e 



On lattice 



[dA^d^d^J ' d3x 'ti 'i(A^) (2) 



T = 



(3) 



with a = a(/3, m) the lattice spacing, = the usual 
lattice variable. The integral in Eq(l) is computed by 
simulating on a lattice with time extension L t and 
space extension with L s 3> L t . 



Renormalization group arguments give 
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(5) 



the (negative) coefficient of the lowest order term of 
the P function. 

The negative sign means asymptotic freedom. 

A j, is the physical scale of the lattice regularized 
QCD. 

T is given by Eq(3) 



T : 
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High T corresponds to small g 2 (order) , low T to large 
g 2 (disorder) , the opposite of what usually happens 
to ordinary systems in statistical mechanics. 

This _peculiar fact naturally brings us to 
Duality!; 



2. Duality 

Duality is a deep concept in statistical mechanics 
and field theory. It applies to systems in [d + 1) di- 
mensions which can have topologically non trivial ex- 
citations in d dimensions. 

These systems admit two complementary descrip- 
tions. 
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1) A direct description in terms of local fields $>(x), 
with (<£>) the order parameters, in which the topolog- 
ical excitations [i are non local . This description is 
convenient in the weak coupling regime g < 1 

2) A dual description in which the topological ex- 
citations /i are local fields and (/i) the (dis)order pa- 
rameters. In this description the original $ fields are 
non local. The dual coupling being gr> ~ ~ this de- 
scription is convenient at large g (strong coupling). 
Duality maps the strong coupling regime of the direct 
description into the weak coupling of the dual and 
viceversa. 

The prototype system for duality is the 2d Ising 
model in which the field is a dicotomic variable a = ±1 
defined on the sites of a square lattice. The La- 
grangean is the sum on nearest neighbors of a fer- 
romagnetic coupling S = —T^ij^aiOj . Putting <? = 7 
the partition function is Z[g, a] — Ee s . 

Identifying one of the coordinate axes ,x, with 
space, the other one t with time the topological ex- 
citations (i are kinks and anti-kinks : a kink fj.(x, t) is 
a spatial configuration at time t with a — —1 at he 
sites 1 < s , a; = +1 at a; > i , an anti-kink has op- 
posite signs. It is easily shown that also the operator 
H which creates a kink has eigenvalues ±1. 

Duality for this model is the equality Q 



where /?' is defined by the equation 

1 



sinh(2/3') = 



sinh(2/3) 



(7) 



(8) 



Eq's(7) and(8) summarize what we have said about 
duality. 

Since in QCD the low temperature phase is disor- 
dered, it is natural to look for dual variables \i and 
for their symmetry , which should be responsible for 
confinement. 

A model example is N = 2 SUSY QCD @ in which 
dual excitations are known to be monopoles. 

Also in ordinary QCD monopoles could be the dual 
excitations 0, and dual superconductivity of the 
vacuum the mechanism of confinement. 

Here the word dual means that the role of elec- 
tric and magnetic fields are exchanged with respect 
to ordinary superconductors : monopoles instead of 
electric charges condense and Meissner effect acts on 
electric field instead of magnetic. 

The pictorial idea is that the chromoelectric field 
acting between a qq pair is channeled into an 
Abrikosov flux tube by dual Meissner effect so that the 
energy is proportional to the distance , which means 
confinement. 

In this mechanism the deconfining phase transition 
is an order-disorder transition from a state with (/x) ^ 
to a state with (/x) = 0,i.e. from superconducting to 
normal. 



3. Lattice QCD 

The only practical known way to study the large 
distance behavior of QCD , which is related to con- 
finement, is to simulate numerically the theory on a 
lattice 

The first question is how to detect confinement and 
dcconfincmcnt on the lattice. The question is non triv- 
ial and parallels the same question in experiments. In 
the absence of quarks (the so called quenched theory) 
, i.e. in pure gauge theory this question has a clear 
answer. In the presence of quarks this is not the case 
any more. 

In quenched theory an order parameter can be de- 
fined [? ] which is the vacuum expectation value of 
the parallel transport along the time direction across 
the lattice, the Polyakov line L(x), which is gauge in- 
variant because of periodic boundary conditions. 



L(x) = Tr[Pe 



(9) 



It can be proved that (L) = e~~ with fi q the chem- 
ical potential of a quark. In the confined phase \x q is 
infinite and {L) = 0. This can be seen alternatively 
as follows. 

Let D(x) = (L'{x)L(u)) be the Polyakov loop cor- 
relator. At large distances by cluster property 



D(x) 



Ce" 



(10) 



On the other hand one has for the potential V(x) 

V(x) = -TlnD{x) (11) 
Together with Eq(10) this gives 

V(x) Wx^od ox (12) 

if (L) = and 



V(x) 



const 



(13) 



if (L) fO. 

(L) is an order parameter for confinement, the cen- 
ter of the group Z3 is the corresponding symmetry. 

A standard finite size scaling analysis of its suscep- 
tibility x = J d 3 xD(x) allows to establish that the 
transition is weak first order [TT| . The transition tem- 
perature it T w 270Mev. 

An alternative order parameter is the vacuum ex- 
pectation value of an operator fi which creates a 
monopole (ll f [Ti| . If the mechanism of dual 

superconductivity is at work, (/i) ^ means dual su- 
perconductivity, (n) — normal vacuum. 

If dual superconductivity is the correct mechanism 
of confinement the behavior of this order parameter 
should coincide with that of the Polyakov line. 

In that case one expects {y) ^ for T < T c , (/i) =0 
for T > T c . (/j) is the ratio of two partition functions 
with the same Boltzman factor /3 = ^ . 
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To detect the phase transition one must go to in- 
finite volume. The dependence on the volume L 3 S of 
susceptibilities is determined by the critical indexes, 
which in turn identify the order of the transition and 
its universality class. 

Instead of the order parameter (p,) itself it proves 
convenient to use the related susceptibility 

P = §p ln (p) ■ 

One has, due to the boundary value (p) = 1 at 

(3 = 

(ri=ef> p V (14) 

For T > T c we find HUE! 

p^-\c\L s + c' (15) 

which by use of Eq(14) means in the infinite volume 
limit (p) — > (normal vacuum). 

For T < T c p — > finite limit as L s — > oo, which, 
again by use of Eq(14) implies (fi) ^ (dual super- 
conducting vacuum). 

For T « T c we expect scaling , since the correlation 
length goes large compared to lattice spacing. 

Dimensional analysis gives 

<A*)«i7/(p^) (16) 

with 7 a possible anomalous dimension, and A the cor- 
relation length of the order parameter. Approaching 
T c , when r = (1 — — > 0, A diverges as 

A OCT"" (17) 

so that j « can be neglected. The variable can 
x 

be traded with tLs by use of Eq.(17) so that 

Qi) n LlgfarL:) (18) 

or 

P ^lIf(tlI) (19) 

A best fit to the data gives for gauge group SU(3) 
v = | which corresponds to first order transition, in 
agreement with the analysis done with the Polyakov 
line quoted above. This indicates that dual supercon- 
ductivity can be the mechanism of color confinement. 

Including dynamical quarks explicitly breaks Z-$ 
symmetry and the Polyakov line is not an order pa- 
rameter any more. 

Moreover there is the phenomenon known as string 
breaking: instead of increasing the potential energy 
when pulling apart a static q—q pair the system prefers 
to create pairs in the form of pions, so that there is 
no string tension but there can be confinement. 

The situation is depicted schematically in Fig.l for 
two quark flavors of equal mass m. At large values 




Figure 1: The phase diagram of N f = 2QCD. The 
transition line is defined by the maxima of the specific 
heat and of the susceptibility of the chiral order 
parameter, m is the quark mass, p the baryon chemical 
potential. 



of the mass quarks decouple and the theory goes to 
quenched . There the phase transition is first order 
and well understood. 

At m = there is a chiral phase transition where 
the spontaneously broken chiral symmetry is restored 
: {ipip) is the order parameter. However at non zero 
values of m chiral symmetry is cxplicitely broken. 

The transition line in the figure is defined by the 
maxima of a number of susceptibilities ( the spe- 
cific heat Cy, the susceptibility of the Polyakov line, 
the susceptibility of the chiral order parameter (tptp)) 
which all coincide within statistical errors. 

The region above the line is conventionally called 
" deconfined" , the region below it " confined" . 

The transition can be explored by use of the dual 
superconductivity order parameter (p) |T3 | and the re- 
sult is that indeed vacuum is a dual superconductor 
in the region below the transition line , and goes to 
normal above it. Moreover the transition is consistent 
with first order at least in the region of low masses. 

The order of the transition across the critical line 
in Fig (1) is a fundamental issue in the study of con- 
finement. As noticed in Sect.l a natural explanation 
of observations on the existence of free quarks is that 
the deconfining phase transition is an order-disorder 
transition : a continuous transition (crossover) would 
imply continuity, and hence would require an unnat- 
ural way of explaining the inhibition factors of order 
w 10~ 15 which are observed in nature. 

An analysis of the chiral transition based on renor- 
malization group and 4 — e can be made |l9j| , assuming 
that the relevant critical degrees of freedom are the 
scalars and pseudoscalars. 

The result is that the chiral transition is first order 
for Nf > 3. 

For Nf = 2 two possibilities exist depending on the 
behavior of axial Ua{X) symmetry across the critical 
point . That symmetry is broken by anomaly at T = 
and is expected to be restored at some temperature. 

If m v ' the mass of the singlet pseudoscalar vanishes 
at T c then the chiral transition is first order and the 
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same is at m ^ in a neighbor of the chiral point [2 
m = 0. [3 

If instead m v > = at T c the chiral transition is 
second order in the universality class of 0(4) and a 
crossover at m ^ 0. In this case a tricritical point ex- [5 
ists at non zero value of the baryon chemical potential [6 
[See. Fig.(l) ] which could be observed in heavy ion 
experiments, but has not been observed up to now. 

The order of the chiral transition can be investi- [8 
gated on the lattice by looking at the behavior of the [9 
specific heat and of other susceptibilities as a function [10 
of tha spatial volume. [11 

This analysis requires large amounts of computer 
time on supercomputers . Pioneering work on the [12 
subject was inconclusive. More recently a big effort [13 
has been put on the problem[2(J which is still going 
on. [14 

The motivation for such an effort is the relevance of 
this issue to the understanding of confinement mecha- [15 
nisms. All the existing evidence points to a first order [16 
transition, or to an order disorder nature of the decon- 
fining transition. This also agrees with the observed [17 
behavior of the order parameter /i for dual supercon- 
ductivity 0|- More work is needed , however , to [18 
fully clarify the situation. 

[19 
[20 

References 

[1] Review of Particle Physics:£PJ15 (2000) 



N. Cabibbo, G. Parisi: Phys. Lett.59B (1975) 67 
R.Hagedom:Suppl.NuovoCimento3(l965) 147 
H.A. Kramers, G.H.Wannier:P/i?/s. Rev.60 (1941) 
252 

L.P.Kadanoff,H.Ceva:P/if/s.i?ei;.B3(1971) 3918 
N.Seiberg,E.Witten:A r McLP/iys.B341(1994) 484 
G.'tHooft in High Energy Physics , EPS Confer- 
ence Palermo 1975. A. Zichichi ed. 
S.Mandelstam:P%;.i?ep.C23(1976)245 
K. Wilson: Phys. Rev. D10(1974) 2445 
AM.Polyakov:Phys.Lett.B72 (1978) 477 
F. Karsch, E. Laerman, C. Schmidt : Phys. Lett. 
B520(2001) 41 

A.DiGiacomo:4cia Phys.Polon. B25 (1994) 215 

L.Dcl Debbio, A.Di Giacomo, G.Paffuti, 

P.Pieri:P%i.Lett.B 355 (1995) 225 

A. Di Giacomo, B. Lucini, L. Montesi, G. Paffuti 

:P%;.i?ei;.D61(2000)034503 

A. Di Giacomo, B. Lucini, L. Montesi, G. Paffuti 

:P/M/s.i?e?;.D61(2000)034504 

A.Di Giacomo,G.Paffuti:P/M/s.#ei;D56(1997) 

6816 

M. D'Elia, A. Di Giacomo, B. Lucini, G. Paffuti, 
C. Pica-.Phys. Rev.B71 (2005) 114502 
R.D.Pisarski,F.Wilczek:P/«/s.i?euD29 (1984) 
338 

M.D'Elia,A.DiGiacomo,C.Pica:P/M/s.i?euD72 
(2005) 114510 



IPM-LHP06-20May 



